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Aims. To investigate the stellar and wind properties of a sample of late-type O dwarfs in order to address the weak wind problem. 
Methods. Far-UV to optical spectra of five Galactic O stars were analyzed: HD 216898 (09IV/08.5V), HD 326329 (09V), HD 66788 
(08V/09V), f Oph (09.5Vnn), and HD 216532 (08.5V((n))). We used a grid of TLUSTY models to obtain effective temperatures, 
gravities, rotational velocities, and to identify wind lines. The wind parameters of each object were obtained by using expanding 
atmosphere models from the CMFGEN code. 
(~ l , Results. We found that the spectra of our sample have mainly a photospheric origin. A weak wind signature is seen in C iv 

^ |. /1/11548,1551, from where mass-loss rates consistent with previous CMFGEN results regarding 08-09V stars were obtained 

I ' (~ 10"'" - 10"' Mq yr"'). A discrepancy of roughly 2 orders of magnitude is found between these mass-loss rates and the val- 

O ' ues predicted by theory (Mvink), confirming a breakdown or a steepening of the modified wind momentum-luminosity relation at log 

Li,ILQ < 5.2. We have estimated the carbon abundance for the stars of our sample and concluded that its uncertainty cannot cause 
^ . the weak wind problem. Upper limits on M were established for all objects using lines of different ions, namely, P v /l/ll 1 18,1 128, 

^ ■ Cm /II 176, N V /11239,1243, Si iv /1/11394,1403, and N iv /11718. All the values obtained are also in disagreement with theoretical 

predictions, bringing support to the reality of weak winds. Together with C iv /1/11548,1551, the use of N v /il239,1243 results in the 
lowest mass-loss rates: the upper limits indicate that M must be less than about -1.0 dex Mvink. Regarding the other transitions, the 
^ ■ upper limits obtained still point to low rates: M must be less than about (-0.5 ± 0.2) dex Mvink- We studied the behavior of the Ho- 

• line with different mass-loss rates. For two stars, only models with very low M's provide the best fit to the UV and optical spectra. We 

' explored ways to fit the observed spectra with the predicted mass-loss rates {Mvink)- By using large amounts of X-rays, we verified 

QQ , that few wind emissions takes place, as in weak winds. However, unrealistic X-rays luminosities had to be used (log Lx/Lggi > -3.5). 

The validity of the models used in our analyses is discussed. 

. Key words, stars: winds - stars: atmospheres - stars: massive - stars: fundamental parameters 
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' 1. Introduction: galaxies by transferring momentum, energy and enriched chem- 

J> ical elements to the interstellar medium (Abbott 1982; Freyer et 

• ^ Massive stars of spectral types O and B play an extremely im- 2003) 

portant role in astrophysics. They possess high effective temper- u u u . j- j f j j *u 

^ ,/^, • . , r Although they have been studied for decades, the properties, 

5—1 atures (Yef/ > lOkK), intense ionizing radiation helds and often . . j i ;• t r^n . . n ^ i J. 

rrt ^ ''.u b ... . origin and evolution of OB stars still present several observa- 

. /T strong mass-losses through stellar winds, making their descrip- i j ,u i u n t-u j j r.u • 

., f ^ , , , • tional and theoretical challenges. The dependency of their mass- 

tion considerably difhcult for atmosphere and stellar evolution , . . n- c i n • 

, , rr., , , • • loss rates (M) on the metallicity (Z) for example, as well as their 

models. These stars are known to be progenitors of fascinat- cp ^- ^ . j • j . » / i ■ \ u 

, . , r. , o • s T ■ T-.1 elfective temperatures and wind structure (e.g. clumping) have 

ing obiects such as Red Supergiants (RSGs), Luminous Blue , . i-. . j • i / 

-ir • 11 /T T,^ r ^ ij- r. /,\ r r. s 11 1 bccn continuousW dcbatcd in thc litcrature duHng thc last ycars 

Variables (LBVs), Wolf-Rayet stars (W-Rs), and thus also of ^.^^ 200^ Martins et al. 2002; Bouret et 

some of the most energetic phenomena in the Universe Le 2005 p^j, 2006; Crowther et al. 2006). 
of type II supernovae and some y-ray bursts (Massey 2003; 

Woosley & Bloom 2006). They also heavily affect their host ^^ong several interesting issues cuiTently under discussion 

(for a review see Puis 2008; Hillier 2008), one that have been 

receiving special attention is the so-called weak wind problem, 

Send offprint requests to: wagner.marcolino@oamp.fr ^j^j^j^ ^^^^^ late-type O dwarf stars. From a qualitative 

* Based on observations made with the NASA-CNES-CSA Far point of view, O stars with weak winds present mainly an absorp- 

Ultraviolet Spectroscopic Explorer and by the NASA-ESA-SERC f- ^ »u u ■ t i ■ a 

, f,, . , ^ _ , ^ , ■ , J, tion spectrum, with the exception being a very few weak wind 

internatjonaj ultra vioJet Explorer, and retrieved from the Multimission ,. ^ , , ^ TiT^^nic-^i-T,^ 

Archive at the Space Telescope Science Institute (MAST). Based on li"^^- 1" f^'^t, often only a weak C iv .1.11548,1551 m P-Cygm is 

observations collected with the ELODIE spectrograph on the 1.93-m seen. In contrast, mid- and early O dwarfs can present P-Cygni 

telescope (Observatoire de Haute-Provence, France). Based on observa- profiles in hnes such as C iv /1/11548,1551, N iv /11718, N v 

tions collected with the FEROS instrument on the ESO 2.2 m telescope, /11239,1243 or O v /11371 (see e.g. Snow et al. 1994; Escolano 

program 074.D-0300 and 075.D-0061. et al. 2008). Quantitatively, weak wind stars are defined as hav- 

** The Appendix is only available in electronic format. ing mass-loss rates of less than about 10"^ Mq yr"'. Bouret et 
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al. (2003) were one of the first to suggest such low values after 
analyzing O dwarfs in the H II region NGC 346 in the SMC. 
The spectra of three objects of their sample could only be repro- 
duced by models using mass- loss rates of ~ 10"'" to 10"'^ Mq 
yr"'. Similar results were found by Martins et al. (2004) for 4 
O dwarfs in the compact star formation region N81, also in the 
SMC. Later, weak winds were also found in some Galactic O 
dwarfs, demonstrating that they are not result of an environmen- 
tal (i.e. Z) effect (Martins et al. 2005). Despite these findings, 
as we will discuss later, very low M values are not of general 
consensus. 

There are interesting questions that O stars with very low 
mass loss rates rise. From the stellar evolution point of view, 
it is well known that mass-loss is a fundamental ingredient in 
the models. Very different values for this parameter can alter 
considerably the way the stars evolve. For example, mass-loss 
can change the rotational structure of a star (decreasing surface 
and internal velocities; Q,) due to removal of angular momen- 
tum combined with internal transport mechanisms (Meynet & 
Maeder 2000). Very low M's might thus imply that stars can 
keep high rotational velocities and get closer to break-up veloci- 
ties. Low mass-loss rates can also change the way we understand 
the evolution in the LBV and W-R phases. The total mass lost 
from the Main Sequence prior the LBV phase can be much less 
than currently thought. As a consequence, in order to be con- 
sistent with observed masses of hydrogen deficient W-R stars, 
other intense mass-loss mechanisms must occur (e.g. continuum 
driven giant eruptions in the LBV phase) and/or evolutionary 
time scales (e.g. of the WNL phase) must be changed (see Smith 
& Owocki 2006; van Marie et al. 2008). Given these considera- 
tions, to be sure that weak winds exist is now essential. 

From the radiative wind theory point of view, the low mass- 
loss rates obtained for late type O dwarfs also present challenges. 
While the last, state-of-art theoretical M predictions of Vink et 
al. (2000; 2001 ; Mvink) present a good match for objects with log 
L/Lq greater than about 5.2 (neglecting clumping), for late-type, 
less luminous objects a discrepancy of even a factor of 100 can 
be found (Martins et al. 2005). Late O dwarfs data also suggest a 
breakdown of the modified wind momentum luminosity relation 
(WLR; Kudrizki & Puis 2000) or at least a change in its slope 
(see Fig. 41 of Martins et al. 2005). All these facts constitute 
what is now called the weak wind problem^ 

A possible reason for the discrepancies aforementioned is 
that mass-loss rates determinations based on UV lines could be 
incorrect, as it was suggested by Mokiem et al. (2007). The use 
of only one UV diagnostic line (i.e. C iv /1/11548, 1551), the abun- 
dance of carbon, the effect of X-rays, and the wind ionisation 
structure derived from the models are claimed to be consider- 
able sources of errors for the derivation of M. This means that 
stronger winds could perhaps have been found if appropriate di- 
agnostic tools were used/available. According to Mokiem et al. 
(2007), three objects which are supposed to be in the weak wind 
regime Oph, CygOB2#2, and HD 217086) actually do not 
show very low mass-loss rates (i.e. M's are consistent with the- 
ory) if Ha is used as the diagnostic. 

In order to clarify some of the issues described above and 
get more insight into the weak wind problem, we have ana- 



' Another type of weak wind problem exists, which is the weaker 
wind signatures of some stars compared to others of the same spectral 
type (e.g. 9' Ori C; Walborn & Panek 1984). Throughout this paper, 
we consider "weaker winds" compared to the predicted mass-loss rates 
according to Vink et al. (2000) and to normal O stars of earlier spectral 
types. 



lyzed in detail far-UV, UV and optical high-resolution spectra 
of five Galactic late-type O dwarf stars. Their stellar and wind 
physical parameters were obtained using the codes TLUSTY 
and CMFGEN (Hubeny & Lanz 1995; HilUer & Miller 1998). 
With this work, the number of 08-9V objects studied by means 
of state-of-art atmosphere models is increased considerably, 
since only six were previously analyzed (Repolust et al. 2004; 
Mokiem et al. 2005; Martins et al. 2005). 

The remainder of this paper is organized as follows. In 
Section |2] we describe the way we have selected our targets and 
the observational material used. A description of the atmosphere 
codes and the adopted assumptions are given in Section [3] The 
analysis of each object of our sample is presented in Section|4] 
The derived stellar and wind parameters are presented later in 
Section |5] along with comparisons to previous results and the- 
oretical predictions. In Section |6] we discuss the carbon abun- 
dance and its relation with the mass-loss rate. In Section |7] we 
present estimates for the mass-loss rates of the programme stars 
using other lines besides C iv /l/l 1548, 1551. In Section [8] we 
discuss the consequences of our findings and explore a way to 
have agreement with the observed spectra using the mass-loss 
rates predicted by theory. Section|9] summarizes the main results 
found in our study. 

2. Target Selection and Observations: 

The Galactic stars claimed to have weak winds in the litera- 
ture belong mainly to the 08, 08.5, 09, and 09.5V spectral 
types (hereinafter we refer to them simply as 08-9V). In order 
to carry out our investigation, we initially selected a sample of 
several objects belonging to these classes by using the Galactic 
O Star Catalog (Maiz- Apellaniz et al. 2004; hereinafter the GOS 
Catalog). Various stars were previously observed by our group 
(at optical wavelengths), and we have retrieved the spectra of 
some others from public available archives. After a first exam- 
ination, we neglected all those objects that presented a peculiar 
spectrum. This included for example stars known to be bina- 
ries (having contamination or known to show strong variability 
in the spectrum) or ON8-9V objects. Due to our aim to study 
wind lines, we have drastically decreased the number of objects 
by considering only the ones having both UV and far-UV data. 
Furthermore, we have also neglected the stars in common with 
the work of Martins et al. (2005), since they were analyzed with 
the same atmosphere code utilized here (CMFGEN). In the end, 
we have chosen a subsample of five objects, listed in Table[T] 

2.1. Optical Data: 

For the optical, we used data collected with the MPI 2.2m tele- 
scope, located at the European Southern Observatory (ESO), in 
La Silla, Chile. The high resolution (R=48000) FEROS spectro- 
graph was used and the total wavelength coverage is ~ 4000- 
9000A. The spectra were wavelength calibrated and optimally 
extracted using the available pipeline (for more details see 
Kauferetal. 1999). 

We have also used spectra obtained with the 1 .93m telescope 
at the Observatoire de Haute-Provence, France, using ELODIE. 
This spectrograph has a resolving power R - 42000 in the wave- 
length range ~3900-6800A. The exposure times were set to 
yield a signal-to-noise (S/N) ratio above 100 at about 5000A. 
The data reduction, the order localization, background esti- 
mate/subtraction, and wavelength calibration, were performed 
using the available pipeline (see Baranne et al. 1996). Each or- 
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far-UV 


uv 


Optical 


Star 


Spectral Type" 


V 




Data Set* 


Date 


Data Set 


Date 


Instrument 


Date 


HD 216898 


09IV, 08.5V 


8.04 


F 


-A05 10303 


2000-08-03 


SWP43934 


1992-02-05 


OHP/ELODIE 


2004-08-29 


HD 326329 


09V 


8.76 


F 


- B0250501 


2001-08-09 


SWP48698 


1993-09-21 


ESO/FEROS 


2005-06-25 


HD 66788 


08V 09V 


9.43 


F 


-P1011801 


2000-04-06 


SWP49080 


1993-11-03 


ESO/FEROS 


2004-12-22 


^Oph 


09.5Vnn 


2.56 




C - C002 


1972-08-29' 


SWP36162 


1989-04-29 


OHP/ELODIE 


1998-06-11 


HD 216532 


08.5V((n)) 


8.03 


F 


-AOS 10202 


2000-08-02 


SWP34226 


1988-09-11 


OHP/ELODIE 


2005-11-10 



^ Spectral types come from Lesh (1968), Schild et al. (1969), Garrison (1970), Walbom (1973), and MacConnell & Bidelman (1976). 

Satellite used to obtain spectra: F=FUSE and C=Copemicus. 
t Orbits #115-247 (Snow et al. 1977). 



der was normalized by a polynomial fit to the continuum, speci- 
fied by carefully selected continuum windows. At last, we have 
merged the successive orders to reconstruct the full spectrum. 
The spectrum of C, Oph was obtained from the ELODIE archive 
(see Moultaka et al. 2004 for more details). We have retrieved 
small portions of its pipeline treated spectrum and then normal- 
ized them using polynomial fits. 

2.2. Ultraviolet Data: 

Different sources were used for the UV and far-UV data. We 
have retrieved spectra from the lUE and FUSE satellites using 
the Multimission Archive at STScfl Regarding the lUE data, the 
high resolution SWP mode was preferred (~ 0.2A), but spectra 
in the LWP/LWR region were also used to help in the determi- 
nation of reddening parameters (R and E(B-V)). Regarding the 
FUSE data, we have relied only on the LiF2A channel, which 
comprises the ~ 1087-1182A interval. Spectra of a same object 
were co-added and smoothed with a five point average for bet- 
ter signal-to-noise ratio (S/N) and clarity. Although small, the 
LiF2A region provided two very useful features for our study, 
namely, P v /i/il 118,1128 and C m /II 176 (see Section Oi. 
We have avoided other FUSE channels covering wavelengths 
less than ~ 1087 A due to severe interstellar contamination. 
For one object of our sample, ( Oph, there is no FUSE data 
available. Therefore, we have used far-UV data obtained with 
the Copernicus ("OAO-3") satellite. The spectrum was retrieved 
from the Vizier Servic^l and normalized "by eye". 

3. Models and Assumptions: 

In order to obtain the stellar and wind parameters of the stars 
of our sample we have used the well known atmosphere codes 
TLUSTY (Hubeny & Lanz 1995) and CMFGEN (Hillier & 
Miller 1998). The TLUSTY code adopts a plane-parallel ge- 
ometry, assumes hydrostatic and radiative equilibrium, and a 
non-LTE treatment including line-blanketing is taken into ac- 
count. As such, it is suitable only to model lines that are not 
formed in a stellar wind, i.e., of photospheric origin. On the other 
hand, the radiative transfer and statistical equilibrium equations 
in CMFGEN are solved in a spherically symmetric outflow. The 
effect of line-blanketing is also taken into account, via a super- 
levels formalism (for more details see Hillier & Miller 1998). 
For the outflow we assume a /? velocity law which connects 
smoothly (at depth) with a hydrostatic structure provided by 
TLUSTY. 

^ MAST - http://archive.stsci.edu/ 

^ http://webviz.u-strasbg.fr/viz-bin/VizieR 



3.1. Methodology: 

To start our analysis, we used a new grid of TLUSTY model 
atmospheres based on the OSTAR2002 grid (Lanz & Hubeny 
2003). The new grid uses finer sampling steps in effective tem- 
perature (lOOO/iT) and surface gravity (0.2 dex), updated model 
atoms (for Ne and S ions, C iii, N ii, N iv, O ii, and O iii), and 
additional ions (Mgn, Aln, Aim, Sin, and Fen; see Lanz & 
Hubeny 2007). In order to derive Te/f, we have used mainly He i 
A44-1 1 and He ii /14542 as diagnostic lines, but fits to other tran- 
sitions such as He i /14713 and He ii /14686 were also checked 
for consistency. Typical eiTors for T,,// range from 1 to 2kK. Our 
derivation of log g was based on the fit to the Hy wings. For this 
parameter, the uncertainty varies from 0.1 to 0.2 dex, depend- 
ing on the object. The rotational velocities (vsini) were adopted 
from previous studies (e.g. Penny 1996; Howarth et al. 1997) 
and/or refined/estimated from the fitting process when neces- 
sary. After we have derived the basic photospheric properties 
(i.e. vsini, log g, and T^jj) we have switched to CMFGEN to 
continue our investigation. As in previous studies, we noted a 
good agreement between TLUSTY and CMFGEN (e.g. Bouret 
et al. 2003). In general, conspicuous discrepancies were sorted 
out when we have increased the number of species or atomic 
levels and/or super-levels in CMFGEN. 

Once effective temperatures were obtained, we have com- 
puted the radii by adopting luminosities typical of the 08-09 
dwarfs (see Table 4 of Martins et al. 2002), following the equa- 
tion Ri, - (L*/47rcrr^^^)'^^. As we did not normalize the UV 
spectra, we have used the reddening parameters R and E(B-V) 
(following Cardelli et al. 1988) and the distance as free parame- 
ters to match the continuum in the lUE region. When a good fit 
was achieved, the R, E(B-V), and the distance used were con- 
sidered representatives for the object in question. The values ob- 
tained were generally consistent with the ones estimated in the 
literature (see Sect. |4]i. Whenever a large discrepancy was en- 
countered, we have revised our adopted and hence keep- 
ing Tt,ff constant. A reasonable agreement is found between the 
final models and the observed absolute fluxes in the UV. In some 
cases however, local scalings of the theoretical continuum had to 
be made. 

The two main wind parameters to be obtained are the mass- 
loss rate and the terminal velocity. Ideally, the determination of 
the mass-loss rates should be done by fitting different (sensi- 
tive) P-Cygni and emission lines in the optical (e.g. Ha) and 
in the UV (e.g. C iv .1.11548,1551, N v .11239,1243, Si iv 
.1/11394,1403). In the case of late-type O dwarfs however, most 
of these wind transitions are absent and the best estimator of M is 
C IV /1/11548,1551. Regarding the velocity structure, we have as- 
sumed P - 1 for all our models. Initial tests with /3 ranging from 
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Table 2. Atomic Data Used in the CMFGEN Models. 



Ion 


Basic Models 
# levels # super-levels 


Full Models 
# levels # super levels 


Hi 


30 


30 


30 


30 


He: 


69 


69 


69 


69 


Hen 


30 


30 


30 


30 


Cii 


22 


22 


338 


104 


Cm 


243 


99 


243 


217 


Civ 


64 


64 


64 


64 


Nm 


287 


57 


287 


57 


Niv 


70 


44 


70 


44 


N V 


49 


41 


49 


41 


On 


296 


53 


340 


137 


Oiii 


115 


79 


115 


79 


Oiv 


72 


53 


72 


53 


Ne II 


- 


- 


242 


42 


Mg II 






65 


22 


Al III 






65 


21 


Si III 






45 


25 


Si IV 


33 


22 


33 


22 


P V 


62 


16 


62 


16 


S IV 






142 


51 


S V 






216 


33 


Fe III 






607 


65 


Fe IV 


1000 


64 


1000 


64 


Fe V 


1000 


45 


1000 


45 



0.8 to 1.6 did not improve the quality of the C iv /}/ll548,1551 
fits. The terminal velocity (voo) is a very difficult parameter to ob- 
tain for 08-09 V stars, since there are no clear, saturated P-Cygni 
profiles. We have chosen to first use previous Voo estimates in the 
literature and then change its value from the fits, if needed. When 
such estimates were not available, we have started our models 
assuming Voo values between 1300-1700 km s~^. We conserva- 
tively assume the uncertainty in this parameter to be about 500 
km s"^ Regarding the mass-loss rate, we consider a conserva- 
tive uncertainty of ~ 0.7 dex, following the analysis performed 
by Martins et al. (2005). 

We have assumed a fixed microturbulence ^, of 15 km s ' in 
the co-moving computation of our models. The CMFGEN code 
allows however a description of a depth dependent microturbu- 
lent velocity when passing to the observer's frame, following the 
equation ^(r) = ^„,/„ -i- {^^a.x - ^mm)v(r)lv^. We have fixed ^,„,„ at 
5 km s"' (near the photosphere) and for ^,„ax we have assumed 
a value of O.Ivk,. As it was already reported in previous studies, 
we noted that slightly different values for these parameters do 
not present significant changes in the spectrum. 

We have used two sets of atomic data for CMFGEN in our 
study. The first fits/models to the observed spectra, tests, and the 
exploration of higher mass-loss rates (Section|2ll were done us- 
ing basic models, i.e., models with a reduced number of species 
and levels/super-levels. For our final model for each object, a 
more complete (full) set was used. In total, 4.5GB (2GB) of 
memory space (RAM) was required to compute each full {ba- 
sic) model. The details of the atomic data are shown in Table 
|2l This approach could save a lot of computational time and re- 
sources, without compromising the results. The main difference 
between these two sets of data is the quality of the optical fit. 
Some lines could only be reproduced after the inclusion of ad- 
ditional species/ions (e.g. Mg n /14481 and Si iii /14553) or by 
increasing the number of levels and/or super-levels. In the ul- 
traviolet, only the 1850-2000A region is affected by the heavier 
data set (being better fitted), due to the inclusion of Fe iii. In other 
parts the changes are minor. A solar chemical abundance was 



adopted for all elements (following Grevesse & Sauval 1998). 
As it will be seen later in Section|4] it provides reasonable fits to 
the observations. In this paper, only the amount of carbon is in- 
vestigated in detail (see Section|6]l, given its direct relation with 
the weak wind problem through C iv /1/11548,1551. 

3.2. Clumping and X-rays: 

Currently, the effect of wind density inhomogeneities (i.e. 
clumping) is implemented in CMFGEN through a depth depen- 
dent volume filling factor / - /t» -H (1 - /t»)e''*''^''^'. The pa- 
rameter \'ci regulates the velocity where clumping starts to be 
important. At velocities higher than v^;, / converges quickly to 
foo- With this description, the model is homogeneous near the 
photosphere. Because we have a lack of wind lines in the spectra 
of the stars studied here, we chose to not use clumping in our 
models, i.e., /oo = 1. We stress that the use of clumping usually 
requires a decrease in the mass-loss rate of a (previously homo- 
geneous final) model by a factor of 1 / ^[fZ, in order to fit (again) 
the observed spectra. Thus, the mass-loss in our final models can 
perhaps be overestimated by a factor of about three if ~ 0.1 
in the wind of these objects. In any case, to not use clumping 
means that we might be biased towards stronger winds. 

CMFGEN also allows the inclusion of X-rays in the models 
to simulate the effects of shocks due to wind instabilities. The 
X-rays emissivity is basically controlled by a shock temperature 
plus a filling factor parameter and it is distributed throughout the 
atmosphere (see Hillier & Miller 1998 for details). The effects 
of X-rays in the atmosphere of O dwarfs were discussed in de- 
tail by Martins et al. (2005). These authors convincingly show 
that X-rays are not as important for early-type as they are for 
late-type O stars (see also Macfarlane et al. 1994). For early- 
type objects (denser winds), the ionisation and wind lines such 
as C IV /1/11548,1551 do not present significant changes when 
X-rays are included. Only wind lines belonging to super-ions 
(O VI /1^1032,1038 and to a lesser extent, N v /11239,1243) are 
affected. In contrast, for late-type stars it is observed that the 
ionisation structure and the C iv /1/11548,1551 line changes con- 
siderably. Because the ionisation is shifted to C v-vi, mass-loss 
rates about ten times stronger than in a model without X-rays are 
sometimes required to fit the observed C iv /1/11548,1551 line. 
Due to this fact, as it is the case of non-clumped models, higher 
mass-loss rates are favored when X-rays are included. In this pa- 
per, we use a fixed shock temperature of 3 x IQ^K and the filling 
factor is chosen to keep log L JLboI close to -7.0 (within +0.1). 
We recall that a log L JLboI of about -7.0 is the typical value ob- 
served (see e.g. Sana et al. 2006). For HD 326329 however, we 
used models with a higher ratio, in accord with the observations 

(log Lx/LboI 6.5; Sana et al. 2006). For ^ Oph, we use a log 

Lx/LboI ~ -7.3 (Oskinova 2005; Oskinova et al. 2006). 

4. Spectral Analysis: 

In this section we present the analysis of our sample. A brief in- 
troduction about each object is made, followed by the presenta- 
tion of our CMFGEN model fits to the UV and optical observed 
spectra. The discrepancies found are discussed and a comparison 
to previous results in the hterature is given when necessary. 

4.1. HD 216898 

HD 216898 belongs to the Cepheus OB3 association (Garmany 
& Stencel 1992). So far, no recent atmosphere models (i.e. uni- 
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Fig. 1. Ultraviolet spectra of HD 216898 (green/light gray line) and our final model (black line; log M - -9.35). 





Fig. 2. Optical spectrum of HD 216898 (green/light gray line) and our final model (black/dashed line). 



fled, Une-blanketed) were used to analyze this object. Its mass- 
loss rate was investigated two decades ago by Leitherer (1988), 
who could only establish an upper limit (log M < -6.97) from 
the Ho- profile. 

Our fits to the far-UV and UV spectra are presented in Figure 
[1] The E(B-V) and distance used are 0.7 and l.lkpc, respec- 



tively. Four intense transitions can be noted in the FUSE region: 
P v AA\ 1 18,1 128, Si iv AA\ 122, 28, and C iii ^1 176. The various 
narrow deep absorptions that are not reproduced by our models 
are of interstellar origin. As it can be seen in Figure [T] (bottom), 
our final model also fits correctly the main features of the lUE 
spectrum. In the 1380-1590A interval. Si iv /l/il394,1403 and 
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Fig. 3. Ultraviolet spectra of HD 326329 (green/light gray line) and our final model (black line; log M - -9.22). A model with a 
higher mass-loss rate is shown as a dotted line (log M - -8.55; see text for more details). 



C IV /1/11548,1551 are the most evident transitions. The other 
several absorptions observed we identified to be mainly due to 
Fe IV lines (an iron forest). We note that a TLUSTY model with 
the same photospheric parameters used in CMFGEN present a 
similar fit to the far-UV and UV, with the exception of C iv 
.i.11548,1551. In TLUSTY, this doublet is seen as two well sepa- 
rated absorptions, while the observations show an extended blue 
feature that is only reproduced with CMFGEN, i.e., it is formed 
in the wind. 

The mass-loss needed to fit C iv /1/11548,1551 is log M - 
-9.35. If we take the other physical parameters derived for this 
object (Teff, L*, and M* = gR^/G) and use the recipe of 
Vink et al. (2000), we derive a predicted mass-loss rate of log 
Mvink - -7.22. A discrepancy of more than two orders of mag- 
nitude is observed ! If we use mass-loss rates higher than the 
one in our final model, a deeper blue absorption (with respect 
to the line center) and an intense emission starts to be seen, i.e., 
a more normal P-Cygni profile is developed, contrary to what 
is observed. On the other hand, lower M's slowly approach the 
photospheric prediction by TLUSTY, as expected. 

In Figure |2] we present our fit to the optical spectrum of HD 
216898. The five different spectral regions shown comprise the 
diagnostic lines used to derive the photospheric parameters (see 
Section [3T1 ). The effective temperature could be very well con- 
strained. The He i ^4471 /He ii ^4542 and He i ^47 13 /He ii 
/14686 ratios are well reproduced with a T^ff = (34 ± l)kK. 
Hy has a reasonable fit but the model is somewhat stronger than 
the observed line. The same happens to Ha. Nevertheless, their 
wings are well matched with a log g of 4.0 ± 0. 1 . Details regard- 
ing this last transition are discussed later in Section 17.21 Some 
of the most intense lines in Figure |2] are: He i /14026, H6, Hy, 
He I /i4388. He i ^4471, He ii ^4542, C iii ^4647, and C iii 
/14650 (with small contribution of O ii). He ii ^4686, He i /14713, 
and Ha. These same transitions are easily identified in the other 
objects of our sample - and in other typical 08-9V stars - that 
have vsini less than about 100 km s"' . For rapid rotators such as 
^ Oph and HD 216532, several profiles are broadened/blended 
and in some cases we cannot distinguish individual transitions 
(see below). We highUght that the Ha observed in HD 216898 
is symmetric. Indeed, all model predictions are symmetrical for 
this line. However, as it is shown later, some objects present an 
asymmetric Ho- profile whose origin is not known. 



4.2. HD 326329 

HD 326329 is located in the NGC 6231 cluster, at the nucleus of 
the Sco OBI association. It belongs to a resolved triple system 
(O9V+O9V-I-B0V; GOS Catalog), but it is probably a single star 
itself since it does not show photometric or spectroscopic varia- 
tions (Garcia & Mermilliod 2001). So far, no atmosphere models 
were used to analyze its spectrum. 

In Figure [3] we present our final model to the far-UV and 
UV observed spectra. From the lUE continua, we have derived 
a E(B-V) of 0.44 and a distance of 1 .99 kpc, in agreement with 
the study of the NGC 6231 cluster made by Baume et al. (1999). 
As in the case of HD 216898, a good agreement is found in the 
FUSE specti-al region. The P v /lil 1 18,1 128, Si iv ^^1122,28, 
and C III /ill76 lines are well reproduced. However, the lUE 
spectrum could not be entirely matched: the observed C iv 
/l/l 1548, 155 1 is considerably deeper than in the model. We have 
tried several different tests to fix this discrepancy (e.g. changes 
in the /3 velocity law; lower terminal velocities; increased tur- 
bulence), but they all turned out to be unsuccessful. We note 
that in the study of Martins et al. (2005) no such problem was 
found. The 09 dwarfs analyzed in their study present a C iv 
/1/11548,1551 profile that was well reproduced by the models. 
We also emphasize that although there is only one lUE high res- 
olution spectrum for HD 326329 (shown in Figure O, model 
comparisons to the other few lower resolution lUE data avail- 
able confirm the too deep/broad feature. Interestingly, one more 
star of our sample also present this characteristic, namely, HD 
216532 (see below). 

Despite the problem described above, we can still be con- 
fident about the mass-loss rate chose for HD 326329, log M = 
-9.22, within ±0.7 dex. First, there is a lack of emission in the 
observed C iv /1/11548,1551 P-Cygni profile, suggesting that M 
must be indeed low. When values for M higher than the one in 
our final model are used, a synthetic line with an intense (un- 
observed) emission and discrepant blue absorption appears (see 
Figure|3j dotted line). On the other hand, when lower mass-loss 
rates are used, C iv /1/11548,1551 slowly turn into photospheric 
features, as the ones predicted in TLUSTY. 

Our fit to the optical spectrum of HD 326329 is presented 
in Figure |4] The He i and He ii lines are well reproduced with 
a Teff = (31 + l)kK. The fit to the Hy wings gives a surface 
gravity (log g) of 3.9 + 0.1. Both T^ff and log g were previously 
estimated by Mathys et al. (2002) who found 31700^: and 4.6, 
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Fig. 4. Optical spectrum of HD 326329 (green/light gray line) and our final model (black/dashed line). 
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Fig. 5. Ultraviolet spectra of HD 66788 (green/light gray line) and our final model (black line; log M = -8.92). 



respectively, by using uvby/3 photometry. Although our derived 
temperature is compatible with their value, their higher log g 
is not supported by our analysis. In fact, a log g much higher 
than 4.0 is not expected for O stars (see e.g. Vacca et al. 1996; 
Martins et al. 2002). Although a fairly good agreement is found 
in several parts of the optical spectrum, we could not achieve 



a good fit to Ha. This line is weaker than in the model and is 
also asymmetric. Its blue wing has a "bump" which distorts the 
profile. This is not seen in the other hydrogen lines. We have not 
used Ha to derive any stellar or wind parameters of the objects 
of our sample. However, we later discuss the effects that different 
mass-loss rates have on this line (see Section|T2]). 
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Fig. 6. Optical spectrum of HD 66788 (green/light gray line) and our final model (black/dashed line). 



4.3. HD 66788 

Figure |5] shows our fits to the ultraviolet spectrum of HD 66788. 
The C IV /1/11548,1551 feature is reproduced with a mass-loss 
rate of log M = -8.92. This value is much lower than the one pre- 
dicted, log Mvink - -6.95. The reddening and distance derived 
are E(B-V) - 0.22 and d = 4.8kpc. Although large, this distance 
is consistent with estimates encountered in the literature (Reed 
1993; Kaltcheva & Hilditch 2000). 

In Figure |6] we show our fit to the optical spectrum. The 
features found in HD 216898 and HD 66788 are very similar, 
suggesting similar photospheric parameters (see Section 14.1b . 
Indeed, from our fits to the Hy we derive a log g = 4.0 ±0.1, 
and from the He i AAAll and He ii /14542 lines we obtain a 
Teff = (34 ± l)kK, as in HD 216898. Theu" rotational velocities 
are also practically the same, about 60 km s"^ We note that the 
Ha line in HD66788 is slightly asymmetric. Although the over- 
all agreement is reasonable, its blue wing has a weak absorption 
not predicted by the model. 



4.4. ^ Oph 

( Oph (also known as HD 149757) is a well studied runaway 
star which presents several interesting characteristics. It has 



been known for a long time to present different kinds of spec- 
tral variability, such as discrete absorption components in the 
UV (DACs), emission line episodes, and line profile variations 
(LPVs) (see e.g. Howarth et al. 1984; Reid et al. 1993; Howarth 
et al. 1993; Jankov et al. 2000; Walker et al. 2005). This object 
presents also a very high rotational velocity: vsini values around 
400 or even 500 km s"' were already reported in the literature 
(Walker et al. 1979; Repolust et al. 2004). 

A reason to include ^ Oph in our sample, is that it is consid- 
ered to have a M estimate based on Wa which is in good agree- 
ment with the predictions of the radiative wind theory (Mokiem 
et al. 2005; 2007). Thus, it is interesting to check if a fit from 
the UV to the optical (including Ha) can be achieved with 
CMFGEN using a low M. We stress however, that the very high 
rotational velocity presented by ^ Oph can bring some problems 
to the analysis. First, fast rotation can distort the shape and in- 
duce temperature and gravity variations through the stellar sur- 
face (see e.g. Fremat et al. 2005). Furthermore, it might imply 
that a stellar wind is not spherically symmetric. In such cases, ID 
unified atmosphere models should be considered as an approxi- 
mation. Another difficulty is that several features in the spectrum 
are broadened and sometimes blended. The analysis of diagnos- 
tic lines thus can present a greater difficulty than in the case of 
slower rotators. 
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Fig. 7. Ultraviolet spectrum of ^ Oph (green/light gray line) and our final model (black line; log M - -8.80). A model with a higher 
mass-loss rate is shown as a dotted line (log M - -8.30; see text for more details). 

Table 3. Stellar and Wind Parameters. 



Star 



HD 216898 HD326329 HD 66788 



^Oph 



HD 216532 



Spec, type 
logg 
T.// (kK) 
vsini (km s"') 
log L/Lo 
R* (Rg) 

M* (Me) 
Voo (km s"') 
log M 

log Mvink 



09IV, 08.5V 
4.0 ± 0.1 
34 ± 1 
60 

4.73 ± 0.25 
6.7_,7 

1700 ± 500 
-9.35 ± 0.7 
-7.22 



09V 
3.9 ±0.1 
31 ± 1 
80 

4.74 ± 0.10 

8.0!l;i 

191 
1700 ± 500 
-9.22 ± 0.7 

-7.38 



08-9V 
4.0 ± 0.1 
34 ± 1 
55 

4.96 ± 0.25 
8.7+3? 

2200 ± 500 
-8.92 ± 0.7 
-6.95 



09.5Vnn 

3.6 ±0.2 

32 ±2 

400 

4.86 ±0.10 

9 2+'-^ 
J3+10 

1500 + 500 
-8.80 ± 0.7 
-6.89 



08.5V((n)) 
3.7 ±0.2 
33 ±2 
190 
4.79 ± 0.25 

'■-^-2.0 

12+5* 
1500 ± 500 
-9.22 ± 0.7 

-6.92 



log LJLboi 
log L/c2 
E(B-V) 
R 

distance (pc) 



-7.00 
-8.42 
0.7 
3.5 
1100 



-6.69 
-8.41 
0.44 
3.1 
1990 



-7.06 
-8.19 
0.22 
2.8 
4800 



-7.31 

-8.29 

0.36 

2.9 

146 



-7.00 
-8.36 
0.72 
3.5 
1150 



Our final model and the lUE spectrum of ( Oph are shown 
in Figure [T] A distance of 146pc and an E(B-V)=0.36 was used 
for the continuum fit. Due to the high vsini involved, the pho- 
tospheric iron forest is considerably broadened. Despite weak, 
N V /II 239, 1243 can be easily distinguished in the spectrum, as 
well as the C ni Alllb and the Si iv /1/11394,1403 lines. As it 
is clear, the model is successful to reproduce the main UV fea- 
tures observed, but some discrepancies can be noted. The C iv 
/1/11548,1551 line for example, could not be matched in detail; 
the synthetic line is slightly more intense and narrower than the 
observed one. A model with a mass-loss rate lower than in our 
final model does not solve the problem and has also a drawback 
in the N v /II 239, 1243 fit. Also, a higher mass-loss rate implies 
in a stronger emission and blue absorption of the P-Cygni pro- 
file, in contrast with the observations (see Figure|2l dotted line). 
The rate used, log M = -8.80, should be certain within a fac- 
tor of three (i.e. ~0.5 dex). This value is much lower than the 
one derived in the study of Mokiem et al. (2005). Using the 
FASTWIND code, and based on Ha, these authors derive log 
M(Ha) = -6.85. Our value is also much lower than the one pre- 
dicted according to Vink et al. (2000), log Mvm - -6.89. We 
come back to this question later in Section]?] Some other smaller 
discrepancies can be also seen in the Si iv /1/11394,1403 line fit 
and on the ~ 1550 - 1650A interval. Despite our efforts, they 
could not be sorted out. Due to the very high rotational velocity 



presented by ^ Oph, future studies using 2D atmosphere models 
might be very useful to address some of these problems. 

In Figure [8] we present our model to the optical spectrum. 
As we mentioned earlier, a high rotational velocity convolution 
is necessary to fit the broad features observed (vsini ~ 400 km 
s"' ). Despite the high rotation, the effective temperature and sur- 
face gravity could be derived, but with a higher uncertainty than 
in the other objects of our sample. From the He i and He ii lines, 
we have obtained a T,,// = (32 + 2)kK. From the Hy profile we 
estimate a low surface gravity compared to the other objects: log 
g = 3.6 + 0.2. It should be kept in mind that this value should 
be regarded as an effective gravity, i.e., the gravity attenuated 
by rotation. A "centrifugal correction" can be applied in order 
to obtain the true gravity, according to the equation log g,rue = 
log igeff + vsini^/Ri,) (see details in Repolust et al. 2004). For 
^ Oph, we find a log gtme of ~ 3.8. This correction however, ne- 
glects any distortion in the stellar shape. A treatment allowing 
the oblateness of the star was presented by Howarth & Smith 
(2001) for three fast rotators, including ^ Oph. These authors 
derived a polar (equatorial) gravity of ~4.0 (~3.6), a polar (equa- 
torial) radius of 7.5Ro (9.IR0), and an effective temperature of 
~ 34kK. The radius obtained from our 7^// and is 9.2Ro. 
Therefore, our stellar parameters for f Oph are compatible to 
the equatorial values found by Howarth & Smith (2001). They 
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Fig. 8. Optical spectrum of ^ Oph (green/light gray line) and our final model (black/dashed line). 



also show agreement with other works in the literature (Repolust 
et al. 2004; Mokiem et al. 2005; Villamariz & Herrero 2005). 



Turning our attention to Ho-, we note that this line and a 
few others were observed to present emission episodes (in the 
form of double peaks) on different occasions (see e.g. Niemela 
& Mendez 1974; Ebbets 1981). However, in general. Ha is ob- 
served to be a symmetric absorption, with an equivalent width 
(Wa) of (2.7 ± 0.4)A, and with a central depth of about 0.8 of the 
continuum intensity. These characteristics can be considered to 
represent ( Oph's quiescent spectrum (see Reid et al. 1993 and 
references therein). In our case, from the observed line we mea- 
sure a Wa of (3.0+0. 3)A and the deepest part of the absorption is 
also at ~0.8 of the continuum level, confirming a quiescent-liks, 
i.e., normal spectrum. However, the profile presents narrower 
wings if compared to the one in our final model (see Figure [8] 
bottom). This problem can be also perceived in the model pre- 
sented by Mokiem et al. (2005) using the FASTWIND code (see 
their Fig. 12). A decrease in the vsini to 350 km s"' improves 
the fit to the wings, but the line center gets slightly deeper than 
observed. Such lower vsini has also a negative effect in other 
optical lines such as Hd and Hy. 



4.5. HD 216532 

HD 216532 is a relatively fast rotator which belongs to the 
Cepheus OB3 association. According to Howarth et al. (1997) 
its vsini is ~190 km s"', a value that we confirm from our 
model fits. Leitherer (1988) have estimated its stellar parame- 
ters and derived an upper limit for the mass-loss rate from Ha, 
log M < -7.08. The present study is the first to analyze quanti- 
tatively its far-UV to optical spectra. 

Our fit to the far-UV and UV data is presented in Figure |9] 
For the E(B-V) and the distance, we inferred 0.72 and 1.15kpc, 
respectively. In the FUSE region, although the overall agree- 
ment is fair, some discrepancies can be noted. The observed 
P V /1.11 118,1 128 and Si iv .1.11122,28 lines are deeper than 
in the model. Although we can see that the synthetic contin- 
uum is somewhat higher than the observed, the problem remains 
if we normalize the spectrum locally. Regarding the lUE spec- 
trum, as in HD 326329, we could not fit the C iv .1.11548,1551 
line. Compared to the profiles found in HD 216898 and HD 
66788, the C iv .1.11548,1551 Une in HD 216532 is broader and 
deeper. Nevertheless, the mass-loss rate of our final model, log 
M = -9.22, should be certain within about +0.4 dex. Higher 
values result in P-Cygni profiles whose emissions are too strong, 
and this is not observed. An example is shown in Figure |9] (see 
the dotted line). As we already mentioned, lower mass-loss rates 
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Fig. 9. Ultraviolet spectra of HD 216532 (green/light gray line) and our final model (black line; log M - -9.22). A model with a 
higher mass-loss rate is shown as a dotted line (log M - -8.78; see text for more details). 



slowly change the C iv /1/11548,1551 feature into two photo- 
spheric absorptions. 

In Figure \W\ we present the observed optical spectrum of 
HD 216532 along with our final model. It can be readily noted 
that the lines are broader if compared to the ones in HD 216898, 
HD 326329, and HD 66788, reflecting the relatively high vsini 
of this object, ~190 km s"'. With an effective temperature of 
~33kK, we achieve a good match to the observed He i AAAll 
and He ii /14542 Unes. To fit the Hy profile, we have used a log 
^ of 3.7. 



5. Analysis of the Results: 

In this section we summarize the results of our spectral analysis 
and make a comparison to previous works and theoretical pre- 
dictions. The stellar and wind parameters obtained for each star 
of our sample are presented in Table[3] For comparison, we also 
list the theoretical mass-loss rates computed following the recipe 
of Vink et al. (2000), Mvink- Additional parameters such as the 
reddening and the distances are shown in the lower part of the 
table. Regarding the X-rays luminosities, we recall that for HD 
216898, HD 66788, and HD 216532, we have fixed log LxILboI 
at the canonical value, i.e. at ~ -7.0. For HD 326329 and ^ 
Oph, we chose LxILboI ratios close to the ones recently observed 
(Oskinova 2005; Oskinova et al. 2006; Sana et al. 2006). 

Overall, the stellar properties of our sample are quite homo- 
geneous. The effective temperatures obtained range from ~30- 
34^7*^, and the radii are between ~1-9Rq. Regarding the surface 
gravities, the most different (lowest) values are presented by HD 
216532 and ( Oph. In both stars, the log g measured should be 
interpreted as effective, i.e., they are attenuated by their fast ro- 
tation. These physical properties show a fair agreement with the 
latest calibration of Galactic O star parameters, regarding the 
08-9V spectral types (see Martins et al. 2005b). 

In order to better analyze the stellar and wind parameters de- 
rived, we also present our results in Figure (TT] in two different 
ways: the mass-loss rates are compared to the theoretical pre- 
dictions of Vink et al. (2000); and the wind parameters are used 
to construct the modified wind momentum luminosity relation 
(WLR): (Mvoc V^) X L*. In each of these plots, we include the 
results of Martins et al. (2005) regarding eleven O dwarfs, as 
well as the data gathered by Mokiem et al. (2007) regarding O 
and B dwarfs, giants, and supergiants (see their Table Al). In or- 



der to simplify the analysis and exclude metallicity effects, only 
Galactic objects are considered. 

In the plots presented in FigurefTT] we can perceive the same 
basic result: the stars of our sample gather around the same 
places occupied by the four 08-9V stars studied in Martins et 
al. (2005), namely, fi Col, AE Aur, HD 46202, and HD 93028. 
This means that: (i) we have also found a discrepancy of roughly 
two orders of magnitude between the measured and the predicted 
mass-loss rates for our programme stars; (ii) the modified wind 
momentum-luminosity relation indeed shows a breakdown or a 
steepening below log L^/Lq ~ 5.2. 

Regarding the brighter objects, i.e. with log L^/L© > 5.2, 
they follow reasonably well the theoretical expectations regard- 
ing the mass-loss rate. We remind that the four most luminous 
early-type dwarfs analyzed in the work of Martins et al. (2005) 
have clumped mass-loss rates. Thus, although they fall some- 
what below the theoretical line in the log M x log Mvink plot and 
below the fit to the data of Mokiem et al. (2007), they do not 
present a significant discrepancy if clumping is neglectecQ. 

The results derived by us and Mokiem et al. (2005) regard- 
ing the star ( Oph are connected in Figure [TT] The large dif- 
ference observed is due to the different mass-loss rates derived. 
While we have log M = -8.80 from the C iv /1/11548,1551 line, 
Mokiem et al. (2005) derived log M - -6.84 based on Hff, us- 
ing the FASTWIND code. The other stellar and wind parameters 
determined in their study and ours are very similar The radius, 
Teff, and log L^^ILq from this (their) paper are 9.2 (8.9) Rq, 32.1 
(32) kK, and 4.86 (4.88), respectively. Regarding the wind ter- 
minal velocity, we (they) obtain 1500 (1550) km s"'. 

As we mentioned above, the objects of our sample and the 
four 08-9V stars analyzed in Martins et al. (2005) occupy about 
the same locus in the panels shown in FigurefTT] It is interesting 
however to analyze their results for other stars having also low 
luminosities, i.e., with log L*/Lq < 5.2. Two objects in Martins 
et al. (2005) have a log Li/L© of about 5.2 (+0.2 dex) and are 
classified as 06.5V stai's: HD 93146 and HD 42088. Another 
object included in their study is HD 152590, an 07.5V with a 
log L*/L0 ~ 4.8. They fall relatively far from the 08-9V group 
in the plots in Figure [TT] Contrary to the other four earlier O 
dwarfs, clumping was not used in their analysis. These three ex- 
amples are suggestive that in stars of the spectral types 06.5V, 



We recall that the theoretical models of Vink et al. (2000) do not 
include clumping. 



12 



Marcolino et al.: Analysis of late-type O dwarfs 




Fig. 10. Optical spectrum of HD 216532 (green/light gray line) and our final model (black/dashed line). 



07V, and 07.5V, a decrease of the wind strength starts to be 
seen, culminating in the 08-9V classes. 



Also at low luminosities are the following objects analyzed 
by Mokiem et al. (2005): Cyg OB2#2 (Bl I), t Sco (B0.2V), 
10 Lac (09V), and HD 217086 (07Vn) (see Figure [TT] filled 
squares). They deserve special attention. In their case, their 
mass-loss rates clearly present better agreement to the theoret- 
ical predictions and also allow them to follow relatively well the 
WLR. Interestingly, their M's were obtained from the Ha line. 
Given these and our UV based results, the origin of the weak 
wind problem can be thought to reside in the different mass-loss 
rates diagnostics employed. In fact, Mokiem et al. (2007) have 
suggested that the problem could be due to uncertainties in the 
UV method. We address and discuss these questions during the 
rest of the paper. 



We emphasize that there are several other 08-9V stars with 
UV spectra similar to the ones found in our sample. Thus, it 
is likely that they possess similar wind properties, i.e., that the 
same basic results we have obtained will be achieved for them if 
we use the same methods of analysis. 



6. The Carbon Abundance 

From the models, it can be easily seen that changes in the carbon 
abundance (ec) affect the C iv /1/11548,1551 profile. In fact, by 
choosing different combinations of M and ec we can get satis- 
factory fits to this feature, as long as its optical depth is kept con- 
stant, i.e., at the observed value. This can be seen more directly 
in terms of the Sobolev optical depth, which is proportional to 
the mass-loss rate, the C iv ionization fraction (qciv), and the 
abundance: tciv M qciv (see for example Lamers et al. 
1999). In the context of the weak wind problem, this suggests 
that we could fit C iv /1/11548,1551 using models with mass-loss 
rates close to the values predicted by the radiative wind theory 
(Vink et al. 2000) and low amounts of carbon. We conclude be- 
low however, that this cannot be the case. 

It is clear that an accurate carbon abundance determination 
in late O dwarfs is needed to disentangle M and ec- In order 
to achieve this goal, we have built a small grid of models with 
the following ec/eco ratios: 0.2, 0.5, 1.0, 1.5, and 3 (mass frac- 
tions). Thereafter, we have analyzed the behavior of the follow- 
ing photospheric transitions: C iii /14647-4651, /II 176, and C iv 
/1/15801,12. Although these lines do respond to changes in ec, an 
accurate determination of the carbon abundance is not feasible 
at the moment. 




Fig. 11. Comparison of our results to previous ones in the literature and theoretical predictions. The objects of our sample are 
represented by star symbols; the O dwarfs studied by Martins et al. (2005) by open triangles; the O and B stars compiled by 
Mokiem et al. (2007) that have log Li,/LQ < 5.2 are represented filled squares and the ones with higher luminosities by filled 
triangles. Left panel: mass-loss rates compared to theoretical predictions (Mvink)- The dashed line indicates a one-to-one relation. 
Right panel: modified wind momentum-luminosity relation (WLR). The solid line represents the best fit to the data in Mokiem et al. 
(2007). The dotted lines link the results obtained in this paper and in Mokiem et al. (2005) regarding the star ^ Oph. Low luminosity 
dwarfs (log L*/Lo < 5.2) in this paper and in Martins et al. (2005) are explicitly indicated by their spectral types (see text for more 
details). 



In Figure [12] we present the fits obtained using models with 
ec = 0.2eQ (20% solai-), ec = (solar), and ec = 3eco (3x 
solar). Although our grid is more complete, we illustrate only 
extreme ec values (besides the solar one) to emphasize the dis- 
crepancies. From Figure[T2] we can note that: (i) when an abun- 
dance of 0.2eca is used, we cannot fit any of the lines shown; (ii) 
when a 3eco value is used we cannot fit properly the C iii /11176 
line, because a large/broad wing appears in the synthetic line and 
this is not observed; (iii) the 3eco abundance is also not favored 
from the fits to C iii ^4647-4651 in HD 216898 and HD 66788, 
and hints to be an excess also in HD326329, ^ Oph, and HD 
216532. From (i) we can conclude that our sample must have a 
ec > 0.2eco- On the other hand, from (ii) and (iii) we are inclined 
to conclude that ec < 3eco- Regarding the C iv /1/15801,12 lines, 
we see that a 3ecQ upper limit is also suggested from the fits. 



From our analysis, we estimate the following abundance 
range for the 08-9V stars in our sample: 0.5 < ec/ec^ < 2. 
Although the uncertainty is large, it is not large enough to allow 
the use of mass-loss rates that are consistent with the values pre- 
dicted by theory (My ink)- A simple test show that if an abundance 
as low as Q-lec^ (not favored by our models) is used along with 
Mvink, the model still present a very intense C iv /l/}1548,1551 
P-Cygni profile which is not seen in any 08-9V stars. Hence, 
although the uncertainty in the carbon abundance adds to the 
uncertainty in M, it cannot cause the weak wind problem. The 
abundance range obtained does not leave enough room for the 
mass-loss rates to be consistent with My ink- A similar conclusion 
was previously achieved by Martins et al. (2005), who quanti- 
fied the error on the mass-loss rate due to the uncertainties in the 
CNO abundances to be about 0.3dex. 



7. Mass-loss Rates in Late-Type O Dwarfs: 

7.1. Other M Diagnostics: 

In most cases, the C iv /1/11548,1551 resonance doublet is the 
only diagnostic available to establish mass-loss rates in late-type 
08-9V stars, which brings uncertainties to the discussion of the 
weak wind problem. In order to overcome this situation, we have 
explored an alternative approach to estimate M values. 

We have proceeded in the following manner: for each ob- 
ject of our sample, we have started with the M determined from 
the best fit to C iv /1/11548,1551 (the final model). Then, we 
have gradually increased this parameter until reach the value pre- 
dicted by the radiative wind theory - Myink - which for each star 
is a function of Te/f, mass and luminosity. We observed that as 
soon as we have used values higher than in our final models (yet 
considerably lower than Mvink), different wind profiles started 
to appear which are not observed. We have thus used this find- 
ing to establish upper limits on M from different transitions, as 
described below. 

First, in Figure [T3] we illustrate the M diagnostics found 
by following the procedure aforementioned. Besides C iv 
/l/}1548,1551, the lines that show significant changes are: P v 
^^1 118,1128, C III Al 176, N v ^1239,1243, Si iv /1/11394,1403, 
and N iv ^1718. When M is low (i.e. ~ lO-^'-lO"^ Mq yr -'), 
N V /1 1239, 1243 is generally weak (or absent) and the other lines 
are essentially photospheric. However, when M is increased they 
gradually have their profiles modified and some develop a very 
intense P-Cygni profile when Mvink is used (see for instance Si iv 
AAl 394, 1403). 

The upper limits (Mupper) for these M-sensitive lines were 
obtained in different ways. For P v ^^1118,1128, C iii ^1176, 
and Si iv /1/11394,1403, we have observed when the profiles 
started to get filled by wind emission and/or a blue absorption 
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Fig. 12. Determination of the carbon abundance in our programme stars. The observed spectra are indicated by solid line lines. The 
models shown have: 20% solar (0.2^^; dotted line), solar (6^; dashed line), and 3x the solar carbon abundance (3eQ;dashed-dotted 
line). 



started to appear, deviating from the (purely photospheric) ob- 
served ones. For N v /11239,1243, we used the strength of the 
whole P-Cygni profile, which generally is weak or absent in 
the observations. The N iv /11718 line has a different response. 
Instead of getting filled in its center, an asymmetric blueward 
absorption profile is formed. We have then determined Mupper 
when a significant (blue) displacement from the observed fea- 



ture was reached. Such peculiar behavior indicates that this line 
is formed in the very inner parts of the stellar wind. In the study 
done by Hillier et al. (2003), a similar situation is found in the 
homogeneous model for the 07 laf-n star AV 83, specially re- 
garding S V /11502 (see their Fig. 10). 

The results obtained for each object of our sample are sum- 
marized in Table|4] The C iv /1/11548,1551 mass-loss rates indi- 



Marcolino et al.: Analysis of late-type O dwarfs 



15 



5 -c 



4 - 




11Q5 



1205 



1405 



1505 



1605 



1705 
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Table 4. Mass-loss Rates and Upper Limits Derived from Different Transitions (log units). 



Star 


M (C IV) 


M (C III) 


M(N V) 


M (N IV) 


M(P V) 


M (Si IV) 


M (Vink) 


HD 216898 


-9.35 


<-7.47 


<- 


iA6 


<-7.80 


<-7.96 


<-7.59 


-7.22 


HD 326329 


-9.22 


<-7.73 


<- 


1.22 


<-7.49 


<-7J3 


<-7.96 


-7.38 


HD 66788 


-8.92 


<-7.22 


<-i 


^.38 


<-7.52 


<-7.52 


<-7.22 


-6.95 


<roph 


-8.80 


<-7.40 


<-i 


^.15 


<-6.89 


<-7.40 


<-7.70 


-6.89 


HD 216532 


-9.22 


<-7.45 


<A 


^.21 


<-7.57 


<-7.45 


<-7.70 


-6.92 



cated are the ones used in our best fits presented in Section|4] The 
predicted mass-loss rates are also shown for comparison. Below 
we present the models used to construct Table H] and comment 
on each object. 



- HD 216898: In Figure \AJi we show the far-UV and UV 
observed spectra of this object along with our final model, 
the models that determine the upper limits in M from each 
line (Mypper', not necessarily the same in each panel), and the 
model with the Vink's mass-loss rate. We can clearly note 
from Figure |AT| that: (i) the final model presents the best fit 
to the observed spectra; (ii) the model with Mvm presents 
the largest discrepancy; and (iii) the models with the upper 
limits for the mass-loss rate do not present satisfactory fits. 
This is also valid for the other objects discussed below. We 
note that when we have increased M, the N v /1 1239, 1243 
line quickly turned to a P-Cygni profile and remained 
essentially unchanged even when we have reached Mvink- 

- HD 326329: An analogous plot as shown for HD 216898 
is presented in Figure IaTSI for this object. We observed that 
the N IV /11718 transition is not so sensitive to M in this 
case. Therefore, the upper limit for M derived from this line 
is quite close to Myink- For the other lines the situation was 
much clearer, i.e., the upper Umits could be easily obtained. 



HD 66788: In Figure [AJ we show the case of HD 66788. 
In general, the same line trends found in HD 326329 and 
HD 216898 are observed. We had no problems in the 
determination of M upper limits from each line. 

^ Oph: The observed spectra and models are shown in 
Figure lA~4l Regarding N iv /11718, we see that our final 
model does not present an adequate fit to the observed line, 
which is considerably deeper. In fact, none of the explored 
M values could fit this feature. By increasing the mass-loss 
rate however, we noted that the absorption still went blue- 
shifted (as in the other objects). We conservatively assume a 
M upper limit from this line equals to the Vink's mass-loss. 
Another difficult situation is presented by N v /1 1239, 1243. 
First, unexpectedly, the model using Mvmk (i e- with the 
strongest mass-loss rate) has a synthetic P-Cygni emission 
less intense than the one in the upper limit model. A slightly 
stronger emission is however found at ~ 1249 A. In order to 
choose the upper limit model shown in the figure, we have 
focused in the absorption part of the P-Cygni profile. The 
theoretical line is deeper than observed, specially around 
1235 A. We have not considered in the analysis the two (too) 
deep absorptions observed as part of the wind profile (see 
the two marks in Figure IA.4I in this line). These features 
present a large variability in f Oph (mainly the bluest one). 
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that the points represent upper limits. 
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Fig. 15. Sensitivity of Ha to different values of the mass-loss 
rate. The model with the lowest mass-loss has log M - -9.22 
and it corresponds to the most intense absorption line. The model 
with the highest mass-loss rate has log M = -6.0 and it corre- 
sponds to the most intense line in emission. The models in be- 
tween are shown by alternating dashed and full lines for clarity. 



which is interpreted as being due to DACs (Howarth et al. 
1993; see their Figure 2). The situation for P v AAl 118,1128 
and Si iv /1/11394,1403 is somewhat simpler. Although 
our final model does not present a very good match to the 
observed lines, it correctly predicts them in absorption. 
The same does not happen with the upper limit and Vink's 
mass-loss models. 



- HD 216532: In Figure lA3] is shown the models and observed 
spectra of HD 216532. We note that the C iii /II 176, P v 
/1.11 118,1 128, and Si iv .1.11394,1403 lines present a well 
developed P-Cygni profile with Mvink- Regarding the P v 
/l/ll 1 18,1 128 lines, the final model fit is not deep enough. 
Nevertheless, an upper limit for M can be derived without 
difficulty as these lines turn to emissions when M is in- 
creased, and this is not what is observed. 

In Figure [14] we explore the results presented above in a 
graphical way. The mass-loss rates obtained are shown relative 
to the predicted value for each star, i.e., relative to Mvm- For 
C IV /1/11548,1551, the mass-loss rates are the ones in our final 
models. For C iii .11176, Si iv .1.11394,1403, N iv .11718, P v 
/l/ll 1 18,1 128, and N v /11239,1243, we remind that the points 
represent only upper limits, i.e., the mass-loss rates must be 
lower than what is displayed. 

Strikingly, we see from Figure[T4]that in all cases the points, 
i.e., upper limits on the mass-loss rate, fall below the theoretical 
line. The only exception is N iv /11718 in the case of ^ Oph, 
for which we have assigned a conservative upper limit value. If 
we consider only the C iv .1.11548,1551 and N v .11239,1243 
lines, we see that most results are below about -1.0 dex, i.e., the 
M's must be less than about ten times the values predicted by 
theory ! If we analyze all other transitions together, we see that 



most points fall at approximately (-0.5 ± 0.2) dex. This means 
that the mass-loss rates are still very low: they must be less than 
about two to five times Mvink ' Taken together, these results bring 
for the first time additional support to the reality of weak winds. 

If we focus only on one object, it is also interesting to note 
in Figure [14] that the points below the highest point (i.e. closest 
to Vink's line) have the respective theoretical lines already in 
disagreement with the observed ones. In this sense, whenever 
we skip to models/points higher than others. Figure [14] can be 
interpreted as errors being accumulated in the model fits. 

We remind that all the models presented in this paper do not 
include clumping. Therefore, some of the upper limits presented 
above could be in fact much lower if the wind of the objects 
studied here are not homogeneous. In this case, we could then 
have that: 



, M clumped , M 
log—r = log- 



Vink 



■ ^og ^/f, 



Vink 



where the first term in the right-hand side is in the vertical 
axis shown in Figure [14] With a typical filling factor / = 0.1, 
some points (depending on the wind clumping sensitivity of the 
respective lines) could go down, i.e. away from the Vink's pre- 
diction, by 0.5 dex. We think however, that it is more reasonable 
to compare the results of homogeneous atmosphere models with 
radiative wind models that are also homogeneous, as it is the 
case in Vink et al. (2000; 2001). 

7.2. Ha behavior: 

Although we have determined mass-loss rates from the far-UV 
and UV spectral regions, we have also analyzed the behavior of 



Marcolino et al.: Analysis of late-type O dwarfs 



17 



Ha to changes in M. This is an important point because as we 
highlighted in Section|5] the use of Hce by Mokiem et al. (2005) 
seems to indicate for a few stars that there are no weak winds, 
i.e., the M values obtained are compatible with Mvink (see also 
the discussion in Mokiem et al. 2007). 

Previously, the work of Martins et al. (2005) had already 
demonstrated two things regarding this line: (i) in a low den- 
sity wind (M ~ 10 '' Mq yr ') the CMFGEN and FASTWIND 
predictions agree; (ii) for M changes from ~ 10"'" to 10"'' Mg 
yr"' the resulting CMFGEN Wa profiles are virtually identical. 
Here, we extend their analysis by presenting the changes in the 
Ha profile for a large range in mass-loss: from ~ 10"'" to 10"'' 
Moyr"'. 

In Figure [TSl we plot models with the photospheric parame- 
ters fixeqj and different values for the mass-loss rate. First, we 
note that only with a log M - -6.30 the Ha line turned to emis- 
sion. For lower rates, the predictions are in absorption and there- 
fore the measurement of M in an observed spectrum must be 
done based in the analysis of the wind filling of the photospheric 
profile. In the example shown, models with a log M = -9.22, 
log M - -7.46, log M - -7.20 do not differ much. Despite a 
factor of one hundred in M, the changes in the profiles are too 
small to allow a secure discrimination of a best fit model to an 
observed spectrum (see below). Overall, the models suggest that 
our confidence in mass-loss rates derived from Ha profiles is 
limited to values > 10"^ M© yr"'. As it can be seen in FigurefTSl 
only beyond this threshold noticeable changes start to be found. 

In Figure [16] we turn our attention to the stars of our sam- 
ple. The fits achieved to Ha using our final models (i.e. with 
M obtained from C iv /1/11548,1551) and models with Mvm are 
presented (see Table[3]for the coiTesponding values). First, if we 
focus on HD 216898 and HD 66788, neither our final models 
nor the ones using Mvwk match the observed Ha. The models 
using Mvink are weaker while our final models are stronger than 
the observed line center. However, the differences compared to 
the observations are very small: ~0.02 of the continuum inten- 
sity. Given the uncertainties involved in the normalization of the 
spectra of our sample around Ha, we find very hard to choose 
with certainty between any of these models (and also the ones in 
between, i.e., with intermediate M's). As in Martins et al. (2005), 
we estimate that the position of the line core can have an error 
up to 2%. However, if a nebular contamination is likely or the 
signal-to-noise ratio of the spectra is not high enough, the uncer- 
tainty increases considerably. Such situations illustrate that it is 
not straightforward to establish M values for 08-9V stars using 
Ha. For HD 326329 the situation is similar, but the Mvmk model 
presents a better fit than our final model. 

Interestingly, for ( Oph and HD 216532, our final models 
seem to be better than the ones using Mvm, and not only in 
the line center. To use the predicted mass-loss rates for these 
stars thus means that we can fit neither the UV nor the opti- 
cal lines! Regarding ^ Oph, Repolust et al. (2004) and Mokiem 
et al. (2005) have presented M measurements from Ha based 
on the FASTWIND code. The former authors derived an upper 
limit of log M - -6.14 and the latter found a mass-loss rate of 
log M = -6.85. In both studies, the parameter used for the 
velocity law was fixed at 0.8. Although the fit in Figure [T6luses 
/3 = 1, we have tested a model with the same predicted mass- 
loss (log Mvink - -6.89) and /3 - 0.8. The synthetic Ha Une 
gets deeper, but it is still not enough to fit the observations. We 
therefore do not confirm the findings of Mokiem et al. (2005). 



^ In this example vsini = 80 km s ' and the others parameters are 
fixed at the values derived for one of our stars, namely, HD 216532. 



We speculate that one reason for this can be the different data 
set used and their reduction. Another alternative is the existence 
of a discrepancy between the predicted Ha profile by CMFGEN 
and FASTWIND at ~ 10"^ Mq yr"'. This is a possibihty since 
in the comparison made by Puis et al. (2005) we can note that 
these two codes can predict slightly different Ha line intensities. 
In the case of an 08V or a OlOV model for example (see their 
Fig. 17), FASTWIND tends to present a stronger Ha absorption 
than CMFGEN. Such differences can affect mass-loss rate mea- 
surements for late O dwarfs obtained from this line. 



8. Discussion: 

It is clear from the previous section that when high mass- 
loss rates are used, i.e. above the values derived from C iv 
/1/11548,1551, a considerable discrepancy is seen between the 
models and the observations (see Figs. IA.lHA.5b . Most impor- 
tantly, this is seen for lines of different ions, for each object of our 
sample. With these different diagnostics, the existence of weak 
winds gains strong additional support. Given the importance of 
this result, we now discuss the validity of our models. 

What are the consequences if we still accept the mass-loss 
rates predicted by Vink et al. (2000) for late O dwarfs as the 
correct ones ? Given our findings, we must conclude that the 
expanding atmosphere models used (from CMFGEN) lack new 
physics or have incorrect assumptions or approximations (or 
both). Even with high mass-loss rates (i.e. equal or close to 
Mvink) the models should provide a reasonable fit to the observed 
spectra, where only a few wind lines are seen. The current wind 
structure/ionisation then must be changed in a drastic way for 
a model to predict only a weak C iv /i/ll548,1551 wind feature 
(and perhaps also a weak N v /11239,1243), despite a high M 
(about one hundred times higher than the ones in this paper or 
even more!). One analogy that can be done is that the uppermost 
synthetic spectrum shown in Figure[T3]would have to turn to the 
one essentially photospheric (bottom). Although possible, this 
scenario presents several difficulties. 

First and foremost, contrary to the recently reported case of 
some B supergiants, where a few key UV wind lines could not be 
reproduced by CMFGEN models (see Searle et al. 2008 and also 
Crowther et al. 20060, we do find agreement with the observa- 
tions (see Figs. [TTfTOli. Although in most spectra only one con- 
spicuous wind signature is observed, the lack of wind lines such 
as P V /l/ll 118,1 128, C iii ^1176, Si iv ^^1394,1403, and N iv 
/11718, is also a constraint attained by our models. In fact, these 
are transitions from the photospheric region (which is smoothly 
connected to the wind) that can be well reproduced. 

In addition, the reasonable fits that we have achieved with 
very low mass-loss rates would have to be merely an unfortu- 
nate coincidence. This would be also true regarding the stars 
analyzed by Martins et al. (2005). Another issue is that depend- 
ing on the physical ingredient(s) neglected or assumptions to be 
revised, the theoretical hydrodynamical predictions will proba- 
bly need to be re-examined. In this case, it would now mean 
that we should trust neither in atmosphere nor in the current ra- 
diative wind models when dealing with late O dwarf stars. The 
weak wind problem would then be ill-defined. Further, any fu- 
ture, more sophisticated models would have to achieve fits to the 



* In this case, a problem in the ionization structure of the models is 
claimed. In fact, their wind structure seems to be distinct from the ones 
found in O stars (Prinja et al. 2005). We note however, that X-rays were 
not used in their analysis and might resolve some of the discrepancies. 
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Fig. 16. Fits to the observed Ha lines in our sample using the mass-loss rates of our final models and the mass-loss rates predicted 
by theory (Mymk)- 



observed spectra similar to ours, with perhaps more free input 
parameters or assumptions (e.g. as in non-spherical winds). 

There is no clear answer as for the missing ingredient(s) 
in the current atmosphere models. For massive stars in gen- 
eral, it is true that although CMFGEN and other codes such as 
FASTWIND are considered state-of-art, different physical phe- 
nomena are not taken into account properly or are not taken into 
account at all due to several technical challenges and lack of ob- 
servational constraints (e.g. non-stationarity; wind rotation; non- 
sphericity; realistic description of clumping; see Hillier 2008 
and references therein). It is possible that advances in some of 
these topics can cast a new light in the analysis of massive stars. 
However, the consequences regarding the weak wind problem 
are difficult to foresee. 

In short, we believe that the inclusion of new physics or re- 
laxation of standard assumptions certainly deserves to be ad- 
dressed in future studies (with the appropriate observational 
constraints), but our atmosphere models for the 08-9V stars are 
the best working hypothesis currently available. 

Below we speculate about some possibilities worth to be in- 
vestigated, and thereafter present the hypothesis of highly (X- 
rays) ionized winds. 

8.1. Magnetic fields and multi-component winds 

From the best fit models, we can examine how closely the mo- 
mentum absorbed by the wind matches that required to satisfy 
the momentum equation. In general, for the 08-9V stars, we 
find that the wind absorbs too much momentum. This can be 
illustrated by noting that the mass-loss rate from a single satu- 
rated line is approximately L/c^ (Lucy & Solomon 1970). In all 



cases, this value exceeds our best fit mass-loss rates (see Table 
[3]l. While none of the wind lines in the models are saturated, 
there are sufficient weak lines to yield a force similar, or larger, 
than the single line limit. 

It is unclear how this discrepancy can be removed. 
Possibilities include the existence of magnetic fields, and/or the 
presence of a significant component of hot gas. The latter could 
be generated by shocks in the stellar wind. As noted by Drew et 
al. (1994) and Martins et al. (2005), it is possible that at the den- 
sities encountered in these winds the shocked gas never cools 
(see also Krticka & Kubat 2009). Such gas would not be eas- 
ily detectable in the UV since it could be collisionally ionized 
to such an extent that very few C IV and N V ions would ex- 
ist. Further evidence for a significant amount of hot gas comes 
from the observed X-ray fluxes; the required filling factors in 
the models (while strongly mass-loss rate dependent) indicate 
that the hot gas in the wind is not a trace component. 

Of interest in this regard are the distinct profiles presented 
by C IV /1/11548,1551. While all show blue shifted absorptions, 
the depth and shapes differ significantly among the stars (see for 
instance the ones in HD 66788, HD 326329, and HD 216532). 
Also of concern is the general weakness of the P-Cygni emission 
components (with the possible exception of ( Oph). Instead of 
very low mass-loss, such an absence might be related to compli- 
cated wind structures. 



8.2. X-rays highly ionized winds: 

Although good fits could be achieved with very low M's, we 
have still studied some alternatives to have agreement with the 
observations using My ink- It became soon clear that the most 
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Fig. 17. Effects of different log Lx I LboI ratios in the ultraviolet synthetic spectrum. Note the similarities between the middle and 
bottom spectra, despite the very different mass-loss rates. 



promising option was to explore high levels of X-rays emission. 
If most of the wind is ionized, i.e. is in the form of a hot plasma, 
we do not expect significant wind emission in the UV lines. 

In order to test the idea above, we explored models with high 
values for the LxILboI ratio. We have found that indeed the wind 
emission is decreased in the desired manner To illustrate this, 
we present in Figure [T7] two models with the same (typically 
Vink) mass-loss rate but with different values for log Lx/LboI- 
the canonical (-7.0) and a much higher ratio (-3.8). In addition, 
a third model with a very low mass-loss rate (log M - -9.3; as 
in a typical final model) is shown with a log Lx/LboI = -7.0. 
Given the high level of X-rays emission, the model with a log 
Lx/LboI - -3.8 has most of its wind ionized and despite the 
high mass-loss rate, no significant line emissions take place, as 
in weak winds ! Indeed, the synthetic spectrum with high M and 
log Lx/LboI resembles the one with very low M and normal log 
Lx/LboI (see the spectra in the middle and bottom of Figure [TtIi. 

In order to better quantify the result above, we have turned 
our attention to the ionisation fractions of P v, N v. Si iv, and 
C IV. For the resonance lines P v AAl 118,1128, N v ^1239,1243, 
Si IV /t/11394,1403, and C iv /1/11548,1551, it is well known that 
what is necessary to fit their observed profiles is an appropri- 
ate value of the product of the mass-loss rate and the respec- 
tive ionisation fraction - M x qi - which means to find their cor- 
rect/observed optical depths (see for instance the case of P v 
AAl 1 18,1 128 in Fullertonet al. 2006). Thus, from our final mod- 
els it can be said that we actually fincQ: 



(1) 



where *P,- is a constant that depends on the observed spectrum 
and qi is the ionisation fraction of the ion /. This latter is usually 
defined as: 



j^'j ni{x)dx 
J^^ n{x)dx 



(2) 



' It is important to note that if the ionisation fractions are reli- 
able/correct, we get M from the model fits. Otherwise, what is actually 
obtained is T;. 



where x - v(r)/vca and n,- and n are the ion and element 
number density, respectively. We can now compute what is the 
ionisation fraction needed to fit the observations by using Mvink, 
which we call qi^vink- Since the constant gives the appropriate 
fit to the observations, we can write: 



(3) 



After we found qi,vi„k from the equation above, for each ion, 
we have computed several models with the mass-loss rate fixed 
at Mvink but with different values of log Lx/LboI- For each model, 
we have derived the ionisation fraction of each ion following 
Eqn. 2. We call these ionisation fractions qx, i.e., the ionisation 
fractions for a specific log Lx/LboI value. We have then verified 
which one of the several qx was equal to qi,vink, in order to find 
for which amount of X-rays we would have agreement with the 
observations using Mvink- 

In Figure [T8]we plot the ratio qiymklqx for each ion versus 
several log Lx/LboI values. For simplicity, we illustrate the case 
of only one object of our sample, HD 216898. The first thing 
to note is that the ions do not behave exactly in the same man- 
ner This is not so surprising, since they have different electronic 
structures and ionization potentials, and thus react differently to 
X-rays. If we focus on C iv and P v, we can see that only when 
very high log Lx/LboI ratios are considered (around -3.0), qx gets 
close to qi.vink- Otherwise, their ratio is much lower than unity. 
In the case of Si iv, qi^vinklqx remains very small regardless the 
value of log Lx/LboI- However, an increase is observed when we 

use a log Lx/LboI near -3.0. For a log Lx/LboI 7 for example, 

we have qi^mklqx ~ 5 x 10"^. For log Lx/LboI 3, this ratio is 

much higher, qiyi„k/qx ~ 0.005. For this ion therefore, the con- 
clusion is that even more X-rays seems to be required to qx reach 
qi.vink- The situation for N v is quite interesting: for both low 
and high values of log Lx/LboI, the qx approaches qiymk- This 
is however, not hard to explain. For high values of log Lx/LboI, 
we observed that practically all nitrogen is concentrated in N vi. 
Thus, the ionisation fraction qx obtained for N v is low as the 
ionisation fraction of Vink, qvink, computed from Eqn. 3. On the 
other hand, when very low log Lx/LboI ratios are considered, we 
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Fig. 18. lonisation fractions as a function of log Lx/LboI- The 
values for qx were computed using Eqn. 2 and Mvink- The qi^vM 
values were derived from Eqn. 3 (see text for more details). 

verified that most of the nitrogen is in N iii-iv. Thus, again, a 
very low qx is obtained for N v. In practice, this means that the 
observed, weak N v /II 239, 1243 line (or its absence!), can be 
reproduced by models using Mvink with either very low or large 
amounts of X-rays. 

Our conclusion from both Figures [17] and [18] is that the 
only way to have high mass-loss rates (Mvink) and still agree 
with the observed spectra is by using a log Lx/LboI k, -3.5. 
Unfortunately, this is not supported by X-rays observations, 
where the log Lx/Lgoi values measured are normally about -7.0 
(see e.g. Sana et al. 2006) with a scatter of about +1.0. It is en- 
couraging however, to explore in future studies different scenar- 
ios of X-rays emission in stellar winds. 

9. Conclusions: 

We have analyzed a sample of five late-type Galactic O dwarfs 
by using atmosphere models from the CMFGEN and TLUSTY 
codes (HD 216898, HD 326329, HD 66788, ( Oph, and HD 
216532). Model fits for the far-UV, UV and optical observed 
spectra were presented, from which stellar and wind parameters 
were obtained. The mass-loss rates were obtained first using the 
C IV /1/11548,1551 feature and then we have explored new di- 
agnostic lines. Our main motivation to study 08-9V stars was 
to address the so-called weak wind problem, recently introduced 
in the literature. The main findings of our study are summarized 
below: 

> The stellar parameters obtained for our sample are quite 
homogeneous (see Table [3]l. Surface gravities of about 
3.8 (after correcting for vsini) to 4.0 and T^j/s of 30 to 
34kK were obtained. These values show a good agreement 
with the latest calibrations of Galactic O star parameters, 
regarding the spectral types of the programme stars. 



> By using the C iv /1/11548,1551 line we have derived 
mass-loss rates considerably lower than theoretical pre- 
dictions (Vink et al. 2000). A discrepancy of roughly two 
orders of magnitude is observed (see Figure [TTI ). We thus 
confirm the results of the study of Martins et al. (2005) 
of weak winds among Galactic 08-9 dwarf stars. We also 
confirm a breakdown or a steepening of the modified wind 
momentum luminosity relation for low luminosity objects 
(logU/Lo <5.2). 

> We have investigated the carbon abundance based on a set 
of UV (C III /11176) and optical photospheric lines. We 
estimated that the following range is reasonable for our 
programme stars: 0.5 < ec/^Co ^ 2 (in mass fractions). 
Although the uncertainty is large, it is not large enough to 
be the reason of the weak wind problem. 

> We have explored different ways (besides using C iv 
/1/11548,1551) to determine mass-loss rates for late O 
dwarfs. We have found that P v /l/ll 1 18,1128, C m /II 176, 
N v /11239,1243, Si iv /1/11394,1403, and N iv /11718 lines 
are good diagnostic tools. We have used each of them to 
derive independent mass-loss rate limits, for each star of our 
sample. We have found that together with C iv /1/11548,1551, 
the use of the N v /1 1239, 1243 line implies in the lowest 
M upper limit rates. The results obtained show that the 
M values must be less than about -1.0 dex compared to 
My ink- By considering the other lines, we still find very low 
mass-loss rates. The results obtained show that M must be 
less than about (-0.5 ± 0.2) dex compared to Mvink- They 
bring additional support to the reality of weak winds. 

> Upper mass-loss rate limits derived from N iv /}1718 and 
C HI /II 176 can be a factor of three, or more, lower than 
the theoretical mass-loss rates of Vink. This is crucial in 
confirming a mass-loss discrepancy as N iv /11718 and to 
a lesser extent C iii /II 176 are formed close to the photo- 
sphere, where the uncertainties in the ionization structure 
should be much less than in the (outer) wind. 

> We have analyzed the Ha line and we observed that its 
profile is insensitive to M changes when we are in the 
< IO^^Mq yr"' regime. For some objects of our sample it 
is uncertain to choose between the Ha fit presented by our 
final models and the ones using Mvink- The interpretation 
of the results is hindered by uncertainties in the continuum 
normalization of echelle spectra. Such situation gets even 
more complicated if a contamination by nebular emission 
is likely. For the stars HD 216532 and ( Oph, the fits to 
the observed Ha lines using Mvink were not satisfactory. 
This result shows that even when Ha is used, lower than 
predicted mass-loss rates can be preferred. 

> We have investigated ways to have agreement with the ob- 
served spectra of the 08-9V stars using the mass-loss rates 
predicted by theory (Mvink)- The only mechanism plausible 
that we found is X-rays. By using high values for the log 
Lx/LboI ratio, we could observe that very few wind emis- 
sions takes place, as in models with very low mass-loss rates 
(~ 10""^ - 10"^ Mq yr"'). However, the values needed to 
be used (log Lx/LboI ^ -3.5) are not supported by the ob- 
servations, which usually measure log Lx/LboI values near 
-7.0. 
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Although our analysis was performed with state-of-art at- 
mosphere models, there are a couple of issues which still need 
to be addressed in future studies. For example, some stars of 
our sample (HD 326329 and HD 216532) present a deep C iv 
/1/11548,1551 feature that could not be well reproduced. The 
origin of this extra absorption remains to be investigated (see 
Section [8T1 ). Studies concerning non-spherical winds, specially 
for stars that are fast rotators, are also of great interest. First 
synthetic spectra computed by our group have shown interest- 
ing line profile changes. A deeper analysis is needed to see if 
some of the discrepancies found in the case of ^ Oph can be 
solved. Furthermore, a more realistic treatment or even alterna- 
tive descriptions of the X-rays emission in the atmosphere mod- 
els (e.g. in conformity with the magnetic confinement scenario; 
Ud-Doula & Owocki 2002) could bring valuable informations. 
First steps in this direction have been taken by Zsargd et al. 
(2009). As we have shown in Section 18.21 X-rays are an effi- 
cient mechanism to change the ionisation structure of the stellar 
winds. 

It would be also very useful to analyze a sample of O stars 
having log L^/Lq near ~ 5.2. The disagreement between at- 
mosphere models and theoretical predictions for Galactic stars 
seems to start around this value. The candidate spectral types 
would be 06.5V, 07V, and 07.5V. Also, this same question 
needs to be better studied in metal-poor environments, i.e. in the 
LMC and SMC. We intend to investigate these issues in a future 
paper 

From the point of view of the hydrodynamics, there are also 
important questions that could be addressed. In the work of Vink 
et al. (1999; 2000; 2001), a predicted mass-loss rate is computed 
based on a set of IS A- WIND models (de Koter et al. 1997) plus 
the use of a Monte-Carlo technique to compute the radiative ac- 
celeration (gi)- Although the results obtained by their procedure 
are self-consistent, it should be kept in mind that the IS A- WIND 
atmosphere models do not include the effect of X-rays. As the 
wind ionization can be changed significantly, new radiative ac- 
celerations can be found and perhaps new predicted mass-loss 
rates might be derived. Within this picture, the agreement with 
the Vink et al. predictions for early O dwarfs may be understood, 
since their winds are not seriously affected by X-rays. 

Our results bring new constraints to the weak wind problem. 
Galactic 08-9V stars seem to present very low mass-loss rates, 
as indicated by the C iv /1/11548,1551 feature, and also addi- 
tional spectral diagnostics. The existence of weak winds poses a 
challenge to the current radiative wind models, and possibly has 
important consequences to the massive stellar evolution theory. 
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Appendix A: M upper limits: 

In this Section, we present the model fits used to establish upper 
Umits on the mass-loss rate from different UV lines. 
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